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Abstract 

The evaluation of the minimal evolution time between two distinguishable states of a system is important 
for assessing the maximal speed of quantum computers and communication channels. Lower bounds for this 
miminal time have been proposed for unitary dynamics. Here we derive an attainable lower bound valid 
for general physical processes, which is connected to the quantum Fisher information for time estimation. 
This result is used to delimit the minimal evolution time for typical noisy channels. 



Introduction 

Quantum mechanics imposes fundamental limits to 
the processing speed of any device as well as to the 
communication speed through any channel. Deriva- 
tion of these basic limits usually assumes that such 
devices are noiseless, undergoing unitary evolutions 
[TH37]. The relevant (and often unwanted) influ- 
ence of the environment on processing or information- 
transfering systems is thus frequently ignored. On 
the other hand, this influence, and in particular the 
decoherence speed, plays an essential role in funda- 
mental physics, especially in the understanding of the 
quantum-to-classical transition [38] . Here we unify 
the description of both computation/communication 
speed and decoherence speed in a single framework, 
which deals with the maximal speed of evolution of 
quantum systems. 

Although much work has been done on the sub- 
ject since the first major result by Mandelstam and 
Tamm 1 , scarce contributions [3131 - I4T] undertake 
nonunitary evolutions. In this Letter, we geometri- 
cally derive a general lower bound to evolution times, 
which depends on the quantum Fisher information 
Fc){t) [42 for the problem of time estimation. Be- 
sides recovering, in the proper limits, previous find- 
ings such as the Mandelstam- Tamm bound [I], our 
result allows the study of experimentally more realis- 
tic open systems and the development of a systematic 



approach for tackling such nonunitary evolutions. We 
exemplify the usefulness of this bound by considering 
typical nonunitary quantum channels 



General bound 

Let D{Fb) be any metric on the space of quantum 
states, assumed to depend only on the Bures fidelity 
Fb, 



tr( 



(1) 



We denote by D{t\,t2) a distance, along this metric, 
between the states of the same system at times t\ and 
t%, i.e., a shorthand notation for D{Fb [p(ti), pfa)]}. 
A bound on D(0,t) is obtained by applying to this 
distance the triangle inequality, considering a division 
of the interval (0, r) into infinitesimal pieces, and us- 
ing the relation between the Bures fidelity and the 
quantum Fisher information |42) . 



F B {ti,h) = l-(i 2 -ii) 



+ 0{t 2 -t 1 f. (2) 



One thus gets a general bound on the distance, 
which depends on the quantum Fisher information 
Jq (t) , and is valid for general physical processes (see 



f 



Supplemental Material for detailed derivation), 



l2d 2 D{F B )/dF B 



[dD{F B )/dF B ] 



(3) 

Notice that this bound is invariant by a rescaling 
D' = kD, and yields an implicit lower bound on the 
evolution time r. 

Although Q is valid for any metric D(Fb), we 
show in the Supplemental Material that the metric 
D(Fb) = 2 arccos V F B (arccos being denned on [0, it] 
throughout the article) leads to the greatest lower 
bound on r, 

2 arccos V^s(0,r) < J yj ' F Q {i)dt, (4) 

where the shorthand notation FB{t\,t-i) for 
Fb [p(ti) , pfe)] is used. One should note that 
the l.h.s. is the shortest distance between the two 
states p(0) and p(r), while the r.h.s. is the length of 
the actual path followed by the state of the system, 
under a given dynamics. This implies, as shown in 
Supplemental Material, that this bound is attained 
if and only if the evolution occurs on a geodesic. 

For a unitary evolution, dictated by the operator 
U(t), the quantum Fisher information is given by 
F Q (t) = 4(AH 2 (t))/h 2 [43], where (AH 2 (t)} is the 
variance in the initial state of an operator H(t) de- 
fined as 



i dt 



(5) 



For a time-independent Hamiltonian H, U(t) — 
e -im/h^ and _ Equation @ yields then 

the Mandelstam-Tamm bound. 

For nonunitary evolutions, the quantum Fisher in- 
formation Fq (t) can be calculated by means of a pu- 
rification procedure [44]. For each system of interest 
S, represented by the density operator ps, we assign 
an environment E. We consider the dynamics of ps 
as resulting from a unitary evolution, corresponding 
to the operator Us,E{t), of a pure state of the en- 
larged system S + E. The quantum Fisher informa- 
tion of S + E, which is an upper bound Cg(i) to the 
quantum Fisher information of subsystem S, is given 



by four times the variance of Hs,E{t)/K defined by 
([5]). There are, in fact, infinitely many different evo- 
lutions of S + E corresponding to the same evolution 
of subsystem S, each of those leading to a value of 
Cq (t) . This freedom is integrally expressed by writing 
the purified unitary evolution as UE(t)Us,E(t), where 
UE(t) is any unitary operator acting on E alone. The 
actual value of Fg(t) for each time t coincides with 
the minimal upper bound Cq (t) over all possible evo- 
lutions of the enlarged system, that is, with the min- 
imization of Cq(£) with respect to UE(t)- Since Cq(£) 
is written in terms of a variance that only depends 
on iiE{t) through h B (t), 



? , . H duL(t) . , . 
i dt 



(6) 



the minimization is performed with respect to 
It should be noted that, since we are interested in 
nonunitary evolutions, the value F B = is typically 
not attained, so that the minimum time corresponds 
in this case to a low, albeit finite, value of F B . 

We consider now examples that illustrate the power 
and usefulness of our general result. 

Amplitude-damping channel 

Let S be a two-state system (states {|0), |1)}), and 
E its environment, which is chosen to start on state 
\0)e- The amplitude-damping channel is described 
by the map 



(7a) 



|0)|0> £ - 

\l)\Q) B ^y/P$)\l)\Q) B - 



y/l-P{t)\Q)\V)E, (7b) 



where the state \1)e is orthogonal to |0)e, and the 
time dependence of P(t) reflects the damping dynam- 
ics. We consider here the paradigmatic exponential 
decay, with rate 7, so that Pit) = e~ 7 *. We note that, 
for the above map, the environment can also be con- 
sidered as a qubit. This channel, which corresponds 
to a nonunitary evolution of S, can be described by 
the unitary evolution operator of S + E, 



Us At) = ex P {-iO(t)(v + d- {E) + 



)}, (8) 



2 



where a± and <Tj_ are raising and lowering operators 
acting respectively on the system and environment 
qubits, and 0(f) = arccos y/P(t). 

The bound can be calculated by taking the variance 
(&H 2 (t)) of H(t) obtained from © via © into ©. 
If S is initially in the ground state, the bound yields 
the trivial result Fg = 1, whereas for S initially in 
the excited state, 

Fb > e-~< T , (9) 

yielding a lower bound for the minimum time r > 
(l/j)ln(l/F B ). The right-hand side of © is the 
actual fidelity corresponding to the channel in this 
case. This shows that this bound is saturated for ev- 
ery r and therefore we have already chosen the best 
purification. Furthermore, the fact that the bound 
is saturated implies that the nonunitary evolution 
connecting the inital and final states is through a 
geodesic path in state space (which includes mixed 
states). This reasoning does not depend on the 
specifics of P(t) and it can be shown (Supplemen- 
tal Material) that the bound is always saturated for 
initially excited states and monotonically decreas- 
ing P(t), which is the relevant behavior for decay- 
ing states [45]. Interestingly, a single-mode Jaynes- 
Cummings (nonmonotonic) dynamics also displays 
such behavior. In this case, S is a two-level atom 
and E is the field, restricted to at most one excita- 
tion. The evolution can be described by Q(t) = gt, g 
being the coupling constant, and obeys ©, with con- 
venient choices of the branches of the square roots. 
Saturation is again observed for every r, without need 
for optimization. These results are atypical, as opti- 
mization is usually necessary for obtaining the high- 
est lower bound for r. 



Markovian dephasing 

Let us now consider Markovian dephasing. System S 
is again a single qubit whose nonunitary evolution is 
described by a map that makes use, as before, of an 



ancilla qubit starting in state |0).e, 

|0)|0) B -> e"*"' (v/P(I)|0)|0) £ + ^~W)\0)\l) E ) , 

|1)|0) E (y/p$}\l)\0) E ^—pJT)\l)\l) E ) , 

(10) 

where Hujq is the energy difference between the qubit 
levels, P(t) := (1 + e~ 7 *)/2, and 7 is the phase-decay 
constant. This description, equivalent to the master 
equation approach, points to the following evolution 
operator: 

U S>E {t) = e'^otZ^^ccos^/pt^ZY^^ ( U ) 

where Z and are Pauli operators acting on the 
system and on the environment qubits, respectively. 

Let %s,E(t) be defined by applying © to the 
evolution operator iiE(t)Us,E{t)', we now minimize 
(AHg E (t)) with respect to fiE(t) from ©, taken 
here as the most general form for a 2x2 Hcrmitian 
operator: 

h E {t) = a (t)xW + (i{t)Y {E ^ + 8{t)Z {E \ (12) 

where the identity is omitted for not changing 
(AM% E (t)). The minimum, written in terms of the 
(constant) average of Z, is 

<A#S,*(t)) = (hW a e-^ + 2jW* ) (f - (Z) 2 ) , 

(13) 

so that © becomes, in terms of elliptic integrals of 
the second kind, E(y,k), 

arccos < \A - (^>V^ + ^ ^ 

with r := cj /7- Eq. (fl"4"l) consistently guarantees the 
eigenstates of Z not to evolve. The corresponding 
bound for the fidelity Fb is compared to an exact 
calculation in Fig. [TJ which shows that the bound 
stays close to the exact result up to the first minimum 
of the latter. 
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Figure 1 : Comparison of the evolution of the bound 
on fidelity Fb (lower surface) as a function of the 
dimensionless time 7T to an exact calculation (upper 
surface) for different values of (Z) (different initial 
states), with r =4. 



In the extreme cases ojq = and 7 = 0, (fT3| yields 
simple analytical expressions for the bound on the 
minimal time, which correspond to the exact results 
when (Z) = 0. For the former, 



1 12 arccos vFb 

r > — In sec 



i-(zy 



(15) 



and for the latter, 



r > 



LUq 



1 



{zy 



(16) 



The bound on the minimal time behaves quite dif- 
ferently in these two situations: for loq = the 
limit guarantees that the evolved and initial states 
never become orthogonal, since even for r — > 00, 
arccos \J Fb reaches, at most, 7r/4; whereas for 7 = 
the bound does not lead to a similar restriction for 
any state other than the eigenstates of Z. There is 
a transition between these two regimes, as seen in 
Fig. [2j which displays the r.h.s. in (fT4|) with r — » 00, 
(Z) = 0, as a function of r. There clearly is a region 
T < r cr - lt where the arccos V Fb does not reach ir/2 



Figure 2: Bound 0{r) for arccos V Fb corresponding 
to single-qubit Markovian dephasing, given by the 
right-hand side of (| X4[) for evolution time r — > 00 and 
(Z) = 0. Values under ir/2 indicate regions where 
the bound guarantees that no state turns orthogonal 
to itself. 

for any finite r, hence the evolved and initial states 
never become orthogonal (for any initial state) . This 
exclusion window - calculations give r cr it — 1.3 - 
yields a simple criterion for defining the regimes of 
strong (r > r cr i t ) and weak (r < r cr i t ) dephasing; 
exclusion windows for reaching given finite fidelities 
can be calculated in an analogous way. 

Minimum evolution time and entanglement 

We now investigate the effect of subsystem correla- 
tions on the evolution speed of a compound system. 
We consider the Markovian dephasing of an N-qubit 
system where each qubit interacts only with its own 
environment, as described by (ITU1) . and compare how 
different initial-state correlations (possibly entangle- 
ment) affect the evolution speed. The evolution op- 
erator is 



he 



N N 
S,E = U E{ 



iV 

'(*)!!{' 



-icjQtZi 



(17) 

-i arccos ^ P(~ft) ZiY^ E) \ 



Since fiE{t) defined by ([6]) now belongs to a 2 N x 2 N 
space, its most general form is rather cumbersome for 
large values of N . We resort instead to minimization 
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over a three-parameter family, hinging on the sym- y 
metry of the system: 

N N N 1 

h E (t) = a(t) ]r if > + m J2 yI e) + s(t) J2 zf\ 

i—l i—l i—l 

(18) 0.001 

where a(t), /3(t), and S(t) are optimization variables. 
We get then 



arccos 





where q := (AZ 2 )/(1 - (Z) 2 ), Z = ^Zj/N, and 
the averages are taken in the initial state. We note 
that < q < 1; for a separable state, g < 1/iV 
(equality if symmetrical on the N qubits) . The val- 
ues q = 1 and (Z) = 0, achievable for GHZ states 
[|0 . . . 0) + e^|l . . . l)]/V2, yield a lower bound valid 
for any possible state. For q = 1, the bound on the 
minimum evolution time scales asr~ l/N through- 
out, see Supplemental Material. 

For separable states, on the other hand, the lower 
bound goes from a r ~ 1/ y/N dependence for 
71//V < wo to ~ l/N for jy/N ^> luq, as shown 
in Fig. [3] and in the Supplemental Material. 

This is a striking result, clearly distinct from the 
one corresponding to unitary evolution. It has al- 
ready been seen in the literature [2"3"H2"T] that for uni- 
tary evolution entanglement is a resource that en- 
hances computation time from a r ~ l/y/N scaling 
(separable, slow state) to t ~ 1/7V (entangled, fast 
state). However, for the nonunitary evolution here 
considered, the minimum evolution time for separable 
states, while scaling with 1/y/N for small AT, eventu- 
ally scales as l/N for jy/N ^> ujq- This result, which 
exemplifies the usefulness of our approach, shows 
that, in the presence of decoherence, the l/N be- 
havior may be attained also for separable states, and 
no matter how small is the decoherence rate 7. One 
should note that the lower bounds obtained here by- 
pass solving the actual equations of motion of the 
system, and are valid for any initial state. 

This novel result can be corroborated for the two 
particular initial states above via direct calculations. 



N 



Figure 3: (Color online) Lower bound (solid curves) 
on time for separable, symmetric state with (Z) = 
to reach Fb = 1%, measured in dimensionless units 
7T as a function of number of qubits N, calculated 
numerically from (fT9"|) . Results from exact calcula- 
tions (|2~T1) are plotted for comparison (dashed curves). 
For the black (blue) curves, r = 4, for the light gray 
(green) curves, r — 200; the asymptotes are propor- 
tional to l/N and l/VN, see Supplemental Material. 



For the GHZ state, the exact fidelity at time r is 
1 + e- N ~i T cos 2Nuj q t 



(20) 



and clearly t ~ l/N for any regime. Meanwhile, for 
the separable, symmetric state with each qubit on the 
equator of its Bloch sphere ((Z) = 0), 



Fr, 



1 



-7T 



cos 2^0''" 



(21) 



and a suitable expansion (Supplemental Material) 
readily shows that for jy/N <tii , r ~ 1/s/N, while 
for 7V A /V" 3> wo, r ~ 1/iV. The same transition to fast 
behavior is then observed, corroborating our previous 
result, see Fig. [3] It should be noted that there is ac- 
tually another constraint to the fast regime, given by 
the breakdown of the Markovian approximation for 
extremely short times. 
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Conclusion 

We have developed a geometrical approach that leads 
to an attainable lower bound for the minimal evolu- 
tion time of dynamical systems. This bound applies 
to both unitary and nonunitary processes, and is ob- 
tained by comparing the actual path followed by the 
system in state space with the distance between the 
initial and final states along a geodesic path, defined 
by a metric that is expressed in terms of the Bures 
fidelity. Whenever the evolution between two states 
is along this geodesic, the bound is tight. Further- 
more, it encompasses several special cases discussed 
in the literature, including unitary evolutions and 
mixed initial states. 

The usefulness of this bound is exemplified by con- 
sidering typical nonunitary quantum channels. For 
the amplitude channel, it leads to a tight bound, 
which evidentiates that the evolution between the 
initial and final orthogonal pure states is along a 
geodesic path though mixed states. For a dephas- 
ing channel, it yields very good lower bounds for the 
mininal evolution time between two non-orthogonal 
states. For N-qubit dcphasing, the evolution speed- 
up due to entanglement of its subsystems, previously 
demonstrated for unitary evolution, is shown to hold, 
in the nonunitary case, also for separable states. 

Our general result allows the estimation of the im- 
pact of the environment on the speed of quantum 
computation and information processing. It is also 
relevant for the estimation of thermalization and de- 
coherence times. 

The authors acknowledge the support of the Brazil- 
ian agencies CNPq, CAPES, FAPERJ, and the Na- 
tional Institute of Science and Technology for Quan- 
tum Information. 
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Supplemental material 



A. Geometric derivation of the general bound 



Let D{Fb) be any function dependent only on the Bures fidelity Fb defined in ([T]). We denote by Fsiti,^) 
the fidelity between the states of the same system at times t\ , t 2 , i.e. a shorthand notation for Fb [p(ti), p(t 2 )]; 
D(ti,t 2 ) will analogously be a shorthand notation for D{Fb [p(ti), p(t 2 )]}. Here we consider dynamical 
evolutions of the system and functions D(Fb) so that D(ti,t 2 ) can be considered as a piecewise smooth 
metric on the space of quantum systems, which implies that dD(t, x)/dx is different from zero when x — > t. 
Dividing a given time interval into an arbitrary number of smaller ones, and using the triangle inequality 
for distances, one gets 



D(0,T)<J2D[(i-l)At,iAt], 



(SI) 



where At = r/n. An expansion in the second argument of the i-th term of the sum around (i — I) At yields 

dD[{i - l)At,x 



D[(i- l)At,iAt] 
so that, when n — ¥ oo, one finds 



dx 



At + 0{At 2 ), 



D(0,t) < f 
Jo 



T dD[t,. 



dx 



(It. 



The integrand can be calculated using the chain rule, 



dD(t,x) 



dx 



dD(F B )dF B (t,x) 



dF E 



dx 



(S2) 



(S3) 



(S4) 



Using Eq. ([2]) of the main text, we see that dFB(t,x)/dx tends to zero when x t. Since D(t,x) is a 
piecewise smooth metric as a function of x, dD(t, x) / dx\ x ^ t is finite, and hence dD(Fs) / 'dFs must diverge 
in the x — > t limit. This indeterminacy can be removed with rHopital's rule, 



dD(t,x) 


dF B (t,x) 
dx 




d 2 F B {t,x) 
dx 2 




dx 


1 


d 


1 


dF B {t,x) 




dD{F B )/dF B 


x^>t dFs 


dD(F B )/dF B _ 


dx 



(S5) 



While the numerator is proportional to Fci(t) due to ([2]), the denominator can be written, by multiplying 
and dividing by dD/dFs, as 



d 2 D(F B )\ fdD(F B )dF B (t, 



[dD{F B )/dF B ] 



dF 2 



dF R 



dx 



(S6) 



The first factor in parentheses can be calculated independently of t, x by replacing the limit x — > t by 
Fb — > 1 since the metric D only depends on t,x through Fs(t,x). This first factor is actually proportional 
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to the curvature of the curve D(Fb) in F B = 1. The second factor in parentheses is simply a recurrence of 
dD(t,x)/dx\ x _ ¥t , the term we are calculating. Substituting in (1S5[) and rearranging the terms, one finds 



dD{t,x) 



-1 2 



<9.i: 



dD{F B ) 



dF B 



d 2 D(F E 
dFl 



(S7) 



Taking the square root of the r.h.s. of the above equation into (|S3[) . and rearranging the terms, one obtains 
the general bound ([3]), valid for any piecewise smooth metric D[F B (t, x)] depending only on the Bures fidelity 

F B (t,x). 

We show now, via purification techniques, that the metric D{F B ) = 2a,rccos^/F B leads to the largest 
lower bound on r. This is done by demonstrating that there exists a path between any two given states po 
and p T , along a geodesic, so that 



TQ(t)dt = 2 arccos y/F B (p ,p T ) ■ 



(S8) 



We start from the particularization of the above relation for pure states under unitary evolution, which was 
demonstrated in [13]: 



Wo) 



sj4(AH2(t))dt = 2 arccos y/ F B {\^ }(^ \, |Vv)(<M), 



(S9) 



where \ip T ) are pure states in a projective Hilbert space and H{t) is the Hamiltonian generating the evolution 
from to \ip T ) along a geodesic. 

Notice now that, to any given initial and final states pa and p T of a system S, one can associate respective 
purifications l^o) and (vPt), so that 

\(^ \^r)\ 2 =F B (p Jr). (S10) 



For these specific purifications, one can now design a Hamiltonian H(t) that connects |\l/o) to \^ T ) along a 
geodesic in the enlarged state space corresponding to S + E. This evolution corresponds to an evolution, in 
general non-unitary, of S, described by the state p(t), which connects po to p T . The corresponding quantum 
Fisher information satisfies, according to [43], the inequality A(/S.H 2 (t)) > J~Q(t), where the variance is 
calculated in the purification 1^(0)). Taking this inequality and (|S10j) into (|S9j) leads to 

/ V -FQ(t)dt < 2 arccos y/F B (po,Pr). (Sll) 
J Po 

Since, from (|4]), the l.h.s. of the above equation must be larger or equal to the r.h.s, it follows that the 
equality sign prevails for the chosen evolution. This implies that the evolution of p(t) is also along a geodesic 
in the corresponding state space. 
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B. Amplitude-damping channel: saturation of the bound for general P(t) 

For the evolution operator Us,E{t) of the amplitude-damping channel displayed in ([Sjl, one has 



dt 



(AHl E (t)) = h 2 (a + a- 



(Sll) 
(Sll) 



For S initially in the ground state, the trivially saturated bound F E > 1 is obtained, whereas an initially 
excited state of S yields 



arccos \l F E < 



f 


de(t) 


dt= [ 


d arccos y/ P{t) 


h 


dt 


Jo 


dt 



dt. 



where we have used Q(t) = arccos \J Pit). Hence, for any monotonically decreasing P(t), 



(S12) 



(S13) 



arccos \/Fb < arccos \J P(t), 
a bound which is clearly saturated by the actual fidelity Fb = P(t). 

C. Dephasing channel: derivation of bound 

It is straightforward to show from the definitions of Hs,E(t) and %s,E{t) that 

n S , E (t) = HsMt) + h'(t), (S14) 

where h'(t) := ul E (t)h E (t)U StE (t). The variance to be minimized is (AW| E (t)), which can be cast in the 
form 

(A#i (JS (t)) = (A#§ )B (t)) + (Ah%(t)) + 2tt[(h'(t)fi s ,E(t)) - (h'(t))(H s , E (t))}- (S15) 

Let | V'o) |0) E be the initial state of S + E. The N-qubit dephasing is described by (fP7|) . from which one 
gets 



ve^ 1 — 1 , 



N 



which imply that 
From (jl8]l . 



{H s ,E{t)) = tHj {Z), 



jVfi 2 7 2 /4 

s 27t _ I 



(h'(t)) = S{t)N(2P(t) - 1) - a(t)N(Z)2y/P(fjy/l-P(t), 
where P(t) := (1 + e-^)/2, and 

(Ah'\t)) =N(a 2 {t) + (3 2 {t) + S 2 {tj) + a 2 (t) (iV 2 (Ai 2 ) - n) 4P(t) (1 - P{t)) 

- N8 2 {t) (2P(t) - l) 2 + 2a{t)8{t)N{Z)2yfplj) y/l - P(t) (2P(i) - 1) . 
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(S15) 
(S15) 

(S16) 
(S17) 

(S18) 



From the previous equations, 

2$l[(h>{t)Hs t E(t)) - (h'(t))(HsMt))} = ~2a(t)huj Q N 2 (AZ 2 )2^Pttj^l - P(t) 



NH-fj3(t) (S19) 
Z), 



Ve 2 ^ - 1 

and (AH 2 S E (t)), according to (|ST5l) . is the sum of (jST6l) . ([Sl8| . and (lST9j) . 
One now minimizes over a(i), <5(i) for each t. From <9(A7i| E (t)}/d(3 = 0, one gets 

/3(t) = ^ 7 (Z). (S20) 

Conditions <9(A7i| E {t))/da — and d(A'Hg E (t)) /dS = lead to a linear system of equations, with 
solutions 



Q(t) = l + ( e 2 ^-l)A g ' (S20) 
fio; e 2 ^(Z)iV g 

^ = - l + (e27t-l)jV g ' (S20) 
where q := (AZ 2 )/(1 - (Z) 2 ). By replacing (IS20l) . (IS20l) . and (IS20)) into (ISTHl) . (ISTHl) . and (|5l9]) . one finds 



' Nqie 2 ^ - 1) + 1 e 2u - 1 

and (TT^l) follows due to (U]). The single-qubit result ([TU)) is recovered by taking N = 1 (note that q = 1 in 
this case). 

D. Dephasing channel: TV-dependence of the most general bound 

The most general bound for N qubits under Markovian dephasing, attained with parameters corresponding 
to the GHZ state, is 



arccos \JFb < VN 

Jo 



N(e 2u - 1) + 1 e 2u - 1 



An equally general, albeit slightly larger, upper bound for arccos s/Fb leads to an analytical expression for 
the bound on r. It is found by considering N(e 2lT — 1) 3> 1. One gets then 



y/N \Jr 2 + 1/4 arctan \J e 2 T T - 1 > arccos \/Fb~, (S23) 

which yields, for N 3> 1, 

1 arccos 2 ^gg 

T -iV 7 (r 2 + l/4)- (S24) 
This expression explicitly exhibits a r ~ 1/iV dependence. An alternative estimate, which leads to a better 
approximation of the integral in (IS22[) . is found in the r ^> 1 limit, yielding, to lowest order in 1/N: 

1 arccos y/F^ ( arccos \/Fb\ , . 

T ~N ujo V + 2r J' (S25) 

Fig EH displays the numerically calculated bound, which is compared to the approximation (|S25[) . 
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1 10 100 1000 10 4 10 5 10 6 



1 10 100 1000 10" 10 5 10 6 



Figure SI: Most general lower bound on time for state to reach Fb — 1%, calculated numerically from 
(|S22[) , as a function of N, with r = 4, 20 and 200, respectively (error bars smaller than symbols). The 
straight line, proportional to l/N, obeys (|S25[) . 



E. Dephasing channel: iV-dependence of the bound for separable, symmetric states 



From (JT9J), the bound for separable states can be written in terms of elliptic functions of the second kind 
E(y,k), 



arccos ^ F B < y 1 - (Z) 2 VN^r 2 + \ 
An expansion to lowest order on jt yields r ~ l/N: 



E\-, 

2' ./,,, 2+ i 



— E arcsine 7T , 



— •yr r 



T > 



1 2 arccos 2 \fF~B~ 
* 7(1 -(Z) 2 ) ' 



For r ^> 1, one gets instead, from (|S26j) . 

arccos ^/Fb~ < yj 1 - (Z) 2 ViVr(l - e" 7T ), 
which leads toar~ l/VN dependence, 

1 arccos \JFb 



T > 



N , 



w ^l-{Zf 



(S26) 



(S27) 



(S28) 



(S29) 



An estimate of the transition between (IS27[) and (|S29|) is obtained by equating the two respective contri- 
butions to the bound on arccos y/Fs- This leads to 



'Ntr ^ r 



7 r tr ~ l/(2r 2 ) 
2 arccos y/Fs 



i - (zy 



(S29) 
(S29) 



resulting in a y/Nt T ~ r scaling. 
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F. Dephasing channel: iV-dependence for separable, symmetric state via exact calculations 

The fidelity of the separable, symmetric state of N qubits, each initially on the equator of its Bloch sphere, 
given by (|2ip . can be rewritten and expanded as 

e-~< T cos(2w r) = 2F B 1/N - 1, (S29) 
(l - IT + + 0( 7 V)) (1 - 2w V + OKV)) = 1 - 2^^ + 0{l/N% (S29) 

where we make use of the fact that r vanishes for increasing values of N . The above expression can be 
simplified to 

1T + M_l! T 2 + 0(7 3 r 3^3 r3) = 2 Hl/F B ) + 0{1/N 2y ^ 

The two regimes can be seen in the above equation: for large values of N the first term of the left-hand 
side is dominant, yielding a r ~ 1/N dependence; for smaller values of N, the second term is dominant, 
and r ~ 1/ VN. We can estimate when this transition occurs by finding the value of r that leads to equal 
contributions from both terms of the left-hand side of (|S30[) . which is r tr = 2j/(4uj 2 — r y 2 ). The corresponding- 
value of N is 

N tI = 2ln(l/ F B )(r 2 -1/4), (S31) 
so that for r> 1 the transition happens around a value N proportional to r 2 . 
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